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A FIRST LOOK WITH Chandra AT SGR 1806-20 AFTER THE GIANT FLARE: 
SIGNIFICANT SPECTRAL SOFTENING AND RAPID FLUX DECAY 
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ABSTRACT 

We report on the results of a ~30 ks Chandra pointing of the soft gamma-ray repeater SGR 1806-20, the first 
X-ray observation with high spectral resolution performed after the 2004 December 27 giant flare. The source was 
found in a bursting active phase and with a significantly softer spectrum than that of the latest observations before 
the giant flare. The observed flux in the 2-10keV range was ~ 2.2 x 10 -11 ergs cm" 2 s" 1 , about 20% lower than 
that measured three months before the event. This indicates that, although its giant flare was w 100 times more 
intense than those previously observed in two other soft gamma-ray repeaters, the post flare X-ray flux decay of 
SGR 1806-20 has been much faster. The pulsed fraction was ~3%, a smaller value than that observed before 
the flare. We discuss the different properties of the post- flare evolution of SGR 1806-20 in comparison to those 
of SGR 1900+14 and interpret the results as a strong evidence that a magnetospheric untwisting occurred (or is 
occurring) after the giant flare. 



Subject headings: stars: magnetic fields - 
X-rays: stars 

1. INTRODUCTION 



stars: pulsars: general — pulsar: individual: SGR 1806-20 



Soft gamma-ray repeaters (SGRs) are neutron stars which 
emit short (< Is) and energetic (< 10 42 erg s" 1 ) bursts of soft 
7-rays. The burst repetition time can vary from seconds to 
years (Gogiis, et al. 2001). During the quiescent state (i.e. out- 
side bursts events) these sources are detected as persistent X- 
ray emitters at a luminosity of ~ 10 35 - 10 36 erg s . Occasion- 
ally, SGRs emit much more energetic "giant flares" (~ 10 44 - 
10 45 ergs _1 ); these are rare events reported until recently only 
on two occasions from SGR 0526-66 and SGR 1900+14 (Mazets 
et al. 1979, Hurley et al. 1999). 

Several characteristics of SGRs, including their bursting ac- 
tivity, are explained in the context of the "magnetar" model 
(Duncan & Thompson 1992, Thompson & Duncan 1995). Mag- 
netars are neutron stars the emission of which is powered by the 
decay of an ultra-strong magnetic field (~ 10 14 - 10 15 G). In this 
model the frequent short bursts are associated to small cracks in 
the neutron star crust, while the giant flares are linked to global 
rearrangements of the star magnetosphere. 

On 2004 December 27 SGR 1806-20 emitted an exception- 
ally powerful giant flare, with an initial hard spike lasting 0.2 
s followed by a ^600 s long pulsating tail (Borkowski et al. 
2004, Hurley et al. 2004, Mazets et al. 2004). The prompt emis- 
sion saturated almost all 7-ray detectors, except for those on the 
GEOTAIL spacecraft which provided a reliable measurement of 
the peak intensity (Terasawa et al. 2005). The isotropic lumi- 
nosity above 50 keV was ^6.47 x 10 47 ergs _1 (for a distance 
of 15 kpc), hundreds of times higher than that of the two gi- 
ant flares previously observed from other SGRs. Following this 
event, afterglow emission, similar to that commonly observed 
in gamma-ray bursts, has been observed in the radio band with 
a resolved extended structure (Cameron et al. 2005, Gaensler et 
al. 2005), and possibly also at hard X-ray energy (Mereghetti 



et al. 2005b). The extremely accurate localization (~0.1") ob- 
tained with the radio data made possible the identification of a 
variable infrared counterpart (Kosugi et al. 2005, Israel et al. 
2005). 

Here we report the results of a Chandra Director's Discre- 
tionary Time observation of SGR 1806-20 which provided the 
first X-ray data set with high spectral resolution after the giant 
flare. 

2. OBSERVATION 

Chandra observed SGR 1806-20 for -30 ks with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) instrument on 2005 
February 8. In order to avoid pile-up and perform pulse phase 
resolved spectroscopy, the source was observed in the Contin- 
uous Clocking (CC) mode, which provides a time resolution 
of 2.85 ms and imaging along a single direction. The source 
was positioned in the back-illuminated ACIS-S3 CCD at the 
nominal target position. Standard processing of the data was 
performed by the Chandra X-ray Center to Level 1 and Level 
2 (processing software DS 7.5.0.1). The data were reprocessed 
using the CIAO software (version 3.2) and the Chandra cali- 
bration files (CALDB version 3.0.0). 

Since in the CC mode the events are tagged with the times 
of the frame store, we corrected the times for the variable delay 
due to the spacecraft dithering and telescope flexure, starting 
from Level 1 data, and assuming that all photons were origi- 
nally detected at the target position 1 . We filtered the data to 
exclude events with ASCA grades 1, 5 and 7, hot pixels, bad- 
columns and possible afterglow events (residual charge from 
the interaction of a cosmic ray in the CCD). In the data pro- 
cessing and analysis we always used the specific bad-pixel file 

1 see the Chandra Science Threads at 

http://asc.harvard.edu/ciao/threads/index.html for details 
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of this observation rather than those provided with the standard 
calibration files. After such filtering the observing time was 
29. Iks. 

3. RESULTS 

3.1. Timing 

In order to carry onto timing analysis we extracted the events 
in the 1-10 keV energy range from a region of 5 x 5 pixels 
around the source position and corrected their arrival times to 
the barycenter of the Solar System. We looked for the pres- 
ence of bursts by binning the counts in intervals of 0.2 s and 
searching for excesses above a count threshold corresponding 
to a chance occurrence of 0.1% (taking into account the total 
number of bins). In this way we identified a single burst, last- 
ing about 0.5 s, at 00:38:25 UT of February 9. 

The relatively poor statistics and small pulsed fraction did 
not permit to determine the pulse period independently from the 
Chandra data. We therefore adopted a period of P=7. 560023 s, 
measured during an almost simultaneous /?XTE observation 
(P. Woods, private communication). For this period the Z 2 test 
(Buccheri et al. 1983) gave a significance of 3.5er for a number 
of harmonics m = 3 (or 2.9cr for m = 2). The resulting pulse 
profiles folded in 16 phase bins are shown in Fig. 1 for three 
different energy ranges (1-10, 1-4, and 4-10 keV). The modu- 
lation is rather low and with some evidence for a double peaked 
profile and, possibly, an energy-dependent shape. By fitting 
with two sine functions the pulse profile in the total energy 
range we obtain pulsed fractions values of PFi = 3.0 ± 1.1% 
and PF 2 = 2.6 ± 1 . 1 % 2 , for the fundamental and for the second 
harmonic, respectively. 



w 

-P 

a 

o 

o 



co 
6 



o 



CD 

d 



in 

CO 

d 



CO 

o 



CO 



CO 

d 




Phase 

FIG. 1 . — Folded pulse profiles in the total (l-10keV; top panel), soft (1- 
4keV; middle panel) and hard (4-10 keV; bottom panel) energy ranges. 



3.2. Spectroscopy 

The source spectrum was extracted from a rectangular region 
of 5 x 25 pixels around the source position, and the background 
independently from a source-free region in the same chip. 

Since the CC mode has not yet been calibrated, standard 
threads for spectral analysis are not available, and the Timed 
Exposure mode response matrices (rmf) and ancillary files (arf) 
are generally used. In order to extract the rmf, we first created 
a weighted map with a re-binning by a factor of 8. We then 
used it in the mkacisrmf tool, with an energy grid ranging 
from 0.3 to lOkeV in 5 eV increments. Using this rmf and the 
aspect histogram created with the aspect solution for this ob- 
servation (asphist), we generated the appropriate arf file for 
the source position. Considering that only a few counts were 
detected below 2 keV, due to the high interstellar absorption, 
and to uncertainties in the instrument calibrations at these low 
energies, we restricted our fits to the 2-8 keV energy range. All 
the fits were performed using Xspec (version 1 1.3). 

Equivalently good results were obtained using either a power 
law (photon index r ~ 1 .8) or a thermal bremsstrahlung model 
(kTbrem ~ lOkeV), while single blackbody and neutron star at- 
mosphere models gave unacceptable fits (x 2 ^2.1 in both cases). 
The results of the acceptable fits are summarized in Table 1, 
where we report for comparison also those obtained in Septem- 
ber 2004, before the giant flare, with XMM-Newton (Mereghetti 
et al. 2005a). The best fit power law spectrum was shown in 
Fig. 2. The absorption derived from the Chandra data was con- 
sistent with the pre-flare value. Keeping the absorption fixed at 

2 we define PF as the semi-amplitude of the two sine functions; all errors in the 
text are at 90% confidence level 



the XMM-Newton value yielded a photon index V = 1. 11 ± 
0.05 and a flux of (2.2±0.2)x 10" 11 ergs cm" 2 s" 1 (2-10 keV, 
not corrected for absorption). 

The addition of a blackbody component to the absorbed power- 
law was not required, contrary to the case of the XMM-Newton 
spectrum of September 2004, which had a higher statistics. Even 
by performing the fit of the Chandra data in the wider 0.3-10 
keV range, the inclusion of an additional blackbody component 
does not improve significantly the fit (F-test probability=0.012). 
Note however that the presence of a blackbody with the temper- 
ature and normalization as seen with XMM-Newton is compat- 
ible with the Chandra data (see Table 1). 

We performed pulse phase resolved spectroscopy by extract- 
ing the spectra for three phase intervals, corresponding to the 
rise and the decay of the broad peak, and to the narrow peak 
(see Fig. 1). The resulting spectral parameters were, to within 
the uncertainties, compatible with those of the phase-averaged 
spectrum. This is not surprising considering that the pulsed 
component represents only ^3% of the total emission. 

We do not find any evidence for absorption or emission fea- 
tures in the source averaged and phase-resolved spectra. We 
derived upper limits as a function of the line energy and width 
(cte) by adding Gaussian lines to the continuum model. For the 
phase-averaged spectrum, the 3er upper limit on the equivalent 
width of narrow lines is ~80 eV. The corresponding values for 
broad lines are 1 10 eV and 135 eV (for <j£=0.1 keV and o-£=0.2 
keV, respectively). 

The burst identified in the Chandra data does not contain 
enough counts for a meaningful spectral analysis. 
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3.3. Extended emission and structures 

In order to search for extended X-ray emission or structures 
around the source, we studied the radial profiles of SGR 1806-20 
Since the CC mode has just one-dimensional imaging capabil- 
ity, we first generated an image of the one-dimensional strip in 
the 0.3-8 keV band, and then subtracted an average count rate 
aimed to remove instrumental and cosmic X-ray background 
(Markevitch & Vikhlinin 2001). 

We then produced a one-dimensional surface brightness dis- 
tribution using the same method as for a two-dimensional ra- 
dial profile 1 . Since calibration Point Spread Functions are not 
available for the CC mode, we extracted in the same way the 
radial profile from a CC-mode observation (Obs-ID 4523) of 
RX J0806+1527. This latter source is known to be point-like 
and has both an X-ray flux and spectrum rather similar to those 
of our target (Israel et al. 2003). 

The ratio between the radial profiles of the two sources, did 
not show any evidence for a significant extended emission or 
structures within 30"around our target, with a 3er upper limit in 
flux of < 10" 14 ergs _1 cm" 2 (in the 4-30"range of radii). 



4. DISCUSSION 

This Chandra X-ray observation of SGR 1806-20 is the 
first with an imaging instrument after the December 2004 gi- 
ant flare. It is therefore interesting to compare the results with 
the pre-flare properties of the source, as measured with XMM- 
Newton in September-October 2004 (Mereghetti et al. 2005a). 

The Chandra data clearly indicate that the spectrum soft- 
ened significantly: we obtained a power law with T ~1.8. This 
has to be compared with the pre-flare values T ~1 .2 (with the 
inclusion of the blackbody) or T ~ 1.5-1. 6 (in the single power 
law model, see Table 1). The flux measured with Chandra 
is ^20% lower than the pre-flare value, but still significantly 
higher than the historical flux level of <~1.3xl0 ergs cm" 2 
s" 1 observed before the second half of 2004 3 . Another differ- 
ence with respect to the pre-flare properties is the smaller pulsed 
fraction (which changed from about 10% to 3%) and the pulse 
profile is now double peaked. 

3 Note that before the giant flare the flux was increasing and the flux level was 
larger than its historical average, see Mereghetti et al. 2005a) 



The post- flare evolution of SGR 1806-20 shows both simi- 
larities and differences when compared to that of SGR 1900+14 , 
the only other case in which good spectral X-ray data have been 
collected after a giant flare. In particular, a significant spectral 
softening was observed to accompany the post-giant flare evo- 
lution of SGR 1900+14 also (Woods et al. 1999; 2001). Even 
though the SGR 1806-20 giant flare was two orders of mag- 
nitude more energetic than that of SGR 1900+14 (and of SGR 
0526-66 as well), it was followed by a very rapid decay of the 
X-ray luminosity. We find that the source flux has dropped be- 
low the pre-flare level after about one month, much faster than 
what observed after SGR 1900+14 giant flare. This suggests 
that the post-flare softening, a feature common to both sources, 
is unrelated to the flare energetics and the decay rate of the X- 
ray flux after the flare. 

SGR 1806-20 and SGR 1900+14 exhibit quite a different 
behavior also in the evolution of their timing properties. The 
pulse profile of SGR 1900+14 changed from a complex, multi- 
peaked pattern to a simpler sinusoidal shape. The pre-flare 
pulse shape has not yet been recovered, a possible signature of 
a permanent rearrangement of the star magnetosphere (Woods 
et al 2001). In the case of SGR 1806-20 the pulse profile 
changed from almost sinusoidal to double-peaked. The data 
obtained during the giant flare indicate a multi-peak structure, 
with a time variable and energy-dependent contribution of the 
different peaks (Palmer et al. 2005, Hurley et al. 2005, Mereghetti 
et al. 2005b). This may indicate a different evolution of the 
geometry of the magnetosphere. Moreover, while the pulsed 
fraction in SGR 1900+14 did not change significantly after the 
flare, we found that in SGR 1806-20 it decreased by about a 
factor of three. 

XMM-Newton observations of SGR 1806-20 carried out few 
months before the giant flare have shown an increase of both the 
spectral hardening and the spin-down rate with respect to his- 
torical values (Mereghetti et al. 2005a). In the picture proposed 
by Thompson, Lyutikov & Kulkarni (2002), a twisted internal 
magnetic field stresses the star solid crust, producing a progres- 
sive increase of the twist angle of the external field lines. A 
giant flare is produced when the crust is not able anymore to 
respond (quasi)plastically to the imparted stresses, and finally 
cracks. The crustal fracturing is accompanied by a simplifica- 
tion of the external magnetic field with a (partial) untwisting of 
the magnetosphere. 
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The spectral softening after the 2004 December 27 event ap- 
pears consistent with such a picture. In fact, the situation after 
the flare is somehow opposite to what occurred before the flare, 
when the twist was increasing. The sudden drop of the external 
twist which followed the giant flare results in a decrease of the 
optical depth to resonant cyclotron scattering in the magneto- 
sphere and hence in a steepening of the power-law spectrum. 

The main observational consequences of a magnetospheric 
untwisting, namely a decrease in the X-ray flux, a softening of 
the spectrum and a decrease of the pulsed fraction (Thompson, 
Lyutikov & Kulkarni 2002; Hurley et al. 2005), appear to be 
present in this first post-flare observation. 
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Table 1 

Spectral results for SGR 1806-20 



Model N H r kT R bb * Flux b xld ( d -°- f -) 

10" 22 cm" 2 keV km ergs _1 cm~ 2 

XMM-Newton Pre-flare (Mereghetti et al. 2005a; Obs. C) 

Power law 6.69±0.13 1.51±0.03 - - 2.65 1.37(72) 

Power law + blackbody 6.51^ 1.21 0.79 ^ 1.9^ 2.65 0.93(70) 



Chandra Post-flare (this work) 



Power law 

Power law 
Thermal Bremsstrahlung 
Power law + blackbody 
Power law + blackbody 



7.1±0.4 1.8±0.1 

6.69 fixed 1.77±0.05 
6.4 ±0.4 

7.5! 2 ^ l-78 + ?g 

6.51 fixed 1.46±0.06 



10 - 7 -2 3 

<0.93 <47.1 
0.79 fixed 1.9 fixed 



2.2 0.76(129) 

2.2 0.76(130) 

2.2 0.78 (129) 

2.2 0.77 (127) 

2.2 0.78 (130) 



a For a distance of 15 kpc; errors in the table are given at 90% confidence level. 

b In the 2-10 keV energy band and in units of 10 -11 ; not corrected for the absorption. 



